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Abstract 

We extend the density matrix renormalization group method to exploit Parity, C<i (ro- 
tation by 7i") and electron-hole symmtries of model Hamiltonians. We demonstrate the 
power of this method by obtaining the lowest energy states in all the eight symmetry sub- 
spaces of Hubbard chains with upto 50 sites. The ground-state energy, optical gap and 
spin gap of regular U = 4i and U = 6t Hubbard chains agree very well with exact results. 
This development extends the scope of the DMRG method and allows future applications 
to study of optical properties of low- dimensional conjugated polymeric systems. 

PACS Numbers: 71.10,+x, 71.28.+d and 78.20.-e 
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The Density Matrix Renormalization Group (DMRG) method [1, 2] has proved to be 
an accurate technique for obtaining a few low-lying states of interacting model Hamiltoni- 
ans with short-range interactions in low-dimensions [3-5]. This technique would find wide 
applications in the electronic structure studies of conjugated polymers and donor- acceptor 
systems, if we can obtain excited eigenstates in different symmetry subspaces. However, 
the DMRG techniques reported in the literature for fermions can exploit only the con- 
servation of the z-component of the total Spin (S^ ot ) and total number of particles(iV tot ) 
in the system[6-8]. Using these symmetries one can obtain a few (about ten) low-lying 
states in different M s and iV sectors [2]. This is due to the limitation on the number of 
eigenstates that can be targetted effectively because of the large matrices that one en- 
counters in the DMRG algorithm. Additionally, the Hamiltonians being spin conserving, 
the states with higher total spin repeat in many M s sectors. For example, triplet states 
will be present in the M s = 0, +1 and —1 sectors. Therefore, obtaining the second singlet 
state usually requires targetting the third state in the M s = sector and eliminating 
the M s = intruder triplet state. Besides, with increasing system size, the number of 
intruders also increases. Thus the inclusion of just these symmetries is not sufficient to 
make detailed comparisons between theory and experiments that probe excited states. 

In system with spatial symmetries, the states important from the point of view of 
optical (one-photon) spectroscopies always lie in a spatial symmetry subspace different 
from the subspace in which the ground state is found. For example, in Hubbard chains, 
the ground state is in the A subspace while the dipole allowed excited states are in the B 
subspace [9]. 

Hubbard chains at half-filling also possess electron-hole symmetry. This allows la- 
belling of the states as either '+', corresponding to covalent subspace, or '— ', correspond- 
ing to ionic subspace. The dipole allowed excited states are found in the '— ' subspace 
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while the ground state is in the '+' subspace. The number of higher energy states in the 
'+' subspace intruding below the lowest '— ' subspace state increases as the strength of 
electron correlation U/t increases [10]. Hence, if this symmetry is not incorporated, many 
'+' space states intrude below the lowest '— ' space state rendering the optically allowed 
excited state difficult to access computationally . 

Obtaining the Hamiltonian in fully block-diagonal form, with respect to the above 
symmetries, is also necessary, if we wish to compute the dynamic nonlinear optic coeffi- 
cients using correction vectors [11]. The correction vectors belong to symmetry determined 
subspaces and are obtained by solving a set of inhomogeneous linear algebraic equations 
which are defined by the Hamiltonian matrix. For example, the correction vector <$~\<jj), 
required for computing first order nonlinearities, is defined by 

{H - Eq + = MG) (1) 

where /tj is the dipole displacement operator, E G and \G) are the ground state energy and 
eigenvector, and u is the excitation frequency. The vector <pf\uj) lies in the l B + space, if 
the Hamiltonian has C2, electron- hole and spin symmetries. If the Hamiltonian matrix is 
not block diagonal then for small values of uj the system is near singular making it almost 
impossible to obtain the solution, (ffp(uj). 

In this communication, we briefly outline a new, symmetrized DMRG procedure which 
exploits fully the above symmetries in the context of DMRG. As an example of its utility, 
we present some excitation gaps for long Hubbard chains and compare them with the 
Bethe- ansatz solutions where possible. 

In the implementation of the symmetrized DMRG scheme, we use the group theoretic 
projection operators for projecting the direct product functions onto different symmetry 
subspaces. The symmetries we have incorporated are the electron-hole symmetry, J, end- 
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to-end interchange symmetry of chains, C 2 and parity, P. The latter two symmetries have 
been employed within the DMRG scheme for spin chains. While the full spin symmetry 
classification can be incorporated into the scheme using Clebsch- Gordon coefficients, for 
the present we have only included parity which bifurcates the space of spin eigenvectors 
into even (e) and odd (o) spin spaces. 

The electron-hole symmetry operator interchanges, with a phase, the creation and 
annihilation operators at a site, 

a\ = {-!)% (2) 

Thus, the Fock space of a single site i, under electron-hole symmetry transforms as Jj|0 >= 
| U>, Ji\ T>= T>, Ji\ i>= (-lY\ i>, Ji\ U>= (-1)|0 >. The full electron-hole 

symmetry operator, J, is given by the direct product of the single site operators. 

J = UJi (3) 

i 

The parity operator, at a site i, Pj, flips the electron spin at a site. Thus, Pj|0 >= |0 >, 
Pi\ i>— I T>, P%\ T>= I i>, P%\ Tl>= ( — 1)1 Tl>- The full parity operator for the system 
is also given by the direct product of the single site parity operators. 

i 

The C*2 symmetry interchanges the states of the left and right halves of the system with 
a phase factor. Thus, C2 operating on the direct product state a, a', p! > gives, 

C 2 \fi, a,a',/i' >= (-l^l^W,/! > ; 7 = (n^ + n eT )(n fl > +n a <), (5) 

where, refers to the p th {p' th ) eigenvector of the left (right) density matrix while 

cr(cr') refers to the Fock space state of the new left (right) site as in the standard DMRG 
procedure, and n M , n^, n a and tv are the occupancies in the states \/i >, |// >, \cr > 
and \a' > respectively. 

4 



The operators C 2 , J and P commute amongst each other and the group formed by 
identity and these symmetries consists of eight elements with the remaining four elements, 
resulting from the closure condition. This group being abelian has eight irreducible rep- 
resentations which are labelled e A + , e A~, °A+, °A~, e B + , e B-, °B + , °B~. 

The projection operator for a given irreducible representation, T, is given by 

Pr = X r (R)R (4) 

R 

where i?'s are the symmetry operations, Xr(R) * s ^ ne character of R in T and h is the order 
of the group. The construction of the symmetry adapted direct product states consists 
in sequentially operating on each of the direct product states by the projection operator. 
The linear dependencies of the symmetry adapted combinations that result are eliminated 
by carrying out a Gram-Schmidt orthonormalization. 

The computational procedure involves obtaining the matrix representation of the sym- 
metry operators of the (2n+2) site chain in the direct product basis. The matrix represen- 
tation of both J and P for the new sites in the Fock space is known from their definitions. 
Similarly, the matrix representation of the operators J and P for the left (right) part of 
the system at the first iteration are also known in the basis of the corresponding Fock 
space states. These are then transformed to the density matrix eigenvectors basis. The 
matrix representation of the symmetry operators of the full system in the direct product 
space are obtained as the direct product of the corresponding matrices, 

< /i, a, a', fi'\R 2n +2W, t, t', v' >=< /i\R n \v >< <j\R\\t X a'\Ri\r' X jj!\R n \v' > (5) 

At the next iteration, we require the matrix representation of the symmetry operators in 
the basis of the eigenvectors of the new density matrix. This is achieved by obtaining the 
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matrix R as a direct product of R n and Hi given by 

< p,a\R n+1 \v,T >=< p\R n \v >< <j\R\\t > (6) 
This matrix is renormalized by the transformation 

R n+1 = O+R^O, (7) 

where O is the matrix whose columns are the chosen eigenvectors of the density matrix. 

The matrix representation of Ci is quite straightforward in that each state \p, a, a', pi > 
is mapped into a state \p', a', a, p > with a phase factor. Thus, each row in the matrix of 
C*2 contains only one element and the entire matrix can be represented by a correspondence 
vector. The matrices J and P are also rather sparse and can be stored in sparse form to 
avoid doing arithmetic with zeros. This aspect of the computation turns out to be crucial 
for the implementation of the scheme since the dimensionality of the direct product space 
is usually very large ( about 10 5 ). 

The coefficients of the direct product functions in the symmetrized basis forms a 
matrix S. The Hamiltonian matrix in the direct product basis can be transformed to the 
Hamiltonian in the symmetrized basis by the transformation 

H2n+2 = S^H 2n+ 2S (8) 

The size of the symmetrized Hamiltonian matrix is of the order ~ 2000 x 2000 for a 
cut-off of m — 100. The low-lying eigenstates of this matrix can be obtained by David- 
son's algorithm. 

The symmetry adaptation scheme described above has been implemented both within 
the infinite chain DMRG algorithm and within the finite system algorithm. In the latter, 
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we incorporate the C 2 symmetry only when the left and right parts of the system are 
identical, i.e. at the end of each finite system iteration. Without iterating over the 
density matrices of the fragments ( i.e. within the infinite system algorithm), we find that 
the energy difference between a chain of length N with N + l and N — 1 electrons is equal 
to the Hubbard parameter U to an accuracy of ~ 1CT 3 . After three iterations of the finite 
system algorithm, the accuracy improves to ~ 10~ 5 , for a U value of At. 

As an example of the use of the above technique, we present results of DMRG calcula- 
tions for uniform Hubbard chains at half-filling, for U/t of 4.0 and 6.0, with chain lengths 
of upto 50 sites. We have obtained the lowest energy states in all the eight subspaces ( 
of (7 2 , J and P) , keeping 70 to 150 eigenvectors of the density matrix. The extrapolated 
ground state energy per site of the infinite chain for both values of U/t agrees with the 
exact values of 0.5737331t and 0.4200125t respectively to 5 decimal places. 

The energy of the one-photon transition from the ground state (lying in the e A + 
subspace) to the lowest energy state in the e B~ subspace defines the optical gap of the 
chain. In fig.l we show a plot of optical gap vs. inverse chain length for U/t = 4.0 and 
6.0 for different values of the cut-off m. The optical gaps for these values of U/t, from the 
Bethe ansatz solution for the infinite chain, are 1.2867/; and 2.8926t respectively [12]. The 
corresponding extrapolated values from a polynomial fit in powers of \/n from DMRG 
are 1.278* and 2.895t, obtained with a cut-off of m=150 for U/t = 4.0 and m=100 for 
U/t = 6.0. We also find that the DMRG optical gap tends to saturate at shorter chain 
lengths as we decrease m. The fit of the optical gap to l/n is a reasonably good straight 
line although the best fit is to a polynomial in this variable. 

In fig. 2 we plot the spin gap which we define as the energy gap between the lowest 
triplet state and the ground state singlet as a function of l/n. The Bethe ansatz solution 
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yields a vanishing spin gap in the thermodynamic limit of the uniform Hubbard chain. 
Polynomial fits to our DMRG data are consistent with this. 

In fig. 3 we present the energy gaps, relative to the ground state, obtained by solving 
for the lowest energy state in each of the eight subspaces, for Hubbard chains with 40 to 50 
sites, for the two values of U/t. The excitations clearly break up into two bands. The lower 
energy excitations correspond to states of different symmetry in the covalent subspace ( 
i.e. predominantly spin excitations) , while the higher energy excitations correspond to 
the ionic subspace ( i.e. charge excitations) . With increasing U/t the states in the 
covalent subspace are less dispersed and so are the states in the ionic subspace. However, 
with increase in U/t the gap between the two bands increases. This feature is in agreement 
with the basic physics of the Hubbard models. 

To summarize, we have presented a symmetrized DMRG scheme for obtaining excita- 
tion gaps in different symmetry subspaces of a large model system. We have illustrated 
our scheme with applications to a uniform Hubbard model, wherein the known excita- 
tion gaps are very well reproduced. We expect that this scheme will prove invaluable in 
modelling the electronic state properties of polymers. 
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Figure Captions: 



Fig.l : 

Optical Gap as a function of inverse chain length for Hubbard chains with U = AM and 
U = 6.01 m corresponds to the number of density matrix eigenvectors retained in the 
DMRG procedure. Arrows indicate the model exact gaps for infinite chains. 
Fig.2 : 

Spin Gap (defined in the text) as a function of 1/n for Hubbard chains with U = AM and 
U = 6M. m corresponds to the DMRG cut-off. Model exact spin gaps vanish for infinite 
chains. 
Fig.3 : 

Energy gaps ( measured from the ground state) of the lowest state in each subspace for 
chain length varying from 40 to 50, for two different values of U/t. The level ordering is 

Ee A + j Eo A + j Eo B + j Ee B + j Ee B - j Ee A ~ j Eo B - j Eo A ~ . 
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